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ABSTRACT 

We use a sample of 45 core collapse supernovae detected with the Advanced Camera for Surveys 
on board the Hubble Space Telescope to derive the core collapse supernova rate in the redshift range 
0.1 < z < 1.3. In redshift bins centered on (z) — 0.39, (z) — 0.73, and (z) = 1.11, we find rates 
3.00±J;gf±J-j$, 7.39±J-J|±f;|g 3 and 9-57±l;lg±^;lg, respectively, given in units of yr^Mpc" 3 10" 4 h%. 
The rates have been corrected for host galaxy extinction, including supernovae missed in highly dust- 
enshrouded environments in infrared bright galaxies. The first errors are statistical while the second 
ones are the estimated systematic errors. We perform a detailed discussion of possible sources of 
systematic errors and note that these start to dominate over statistical errors at z > 0.5, emphasizing 
the need to better control the systematic effects. For example, a better understanding of the amount of 
dust extinction in the host galaxies and knowledge of the supernova luminosity function, in particular 
the fraction of faint M ^ — 15 supernovae, is needed to better constrain the rates. When comparing 
our results with the core collapse supernova rate based on the star formation rate, we find a good 
agreement, consistent with the supernova rate following the star formation rate, as expected. 

Subject headings: supernovae: general - surveys 



1. INTRODUCTION 

Supernovae (SNe) mark the end of the life cycle of 
certain stars. Studies of SNe are important both for un- 
derstanding the physics leading to and driving these ex- 
plosions, as well as the impact that the SNe have on their 
environments. As the main producers of heavy elements, 
SNe are pivotal for understanding the chemical evolution 
in galaxies as well as the intergalactic medium via SN- 
driven outflows from galaxies. SNe are also thought to 
be one of the producers of dust in the universe (e.g., re- 
view by Kozasa et al. 2009). Furthermore, SNe tie into 
feedback processes regulating galaxy formation. To un- 
derstand how these processes evolve with cosmic time, it 
is important to understand how the number of explod- 
ing SNe changes with redshift, i.e., the evolution of the 
cosmic SN rate (SNR). Since SNe, particularly the core 
collapse (CC) SNe, have massive progenitors with short 
main-sequence lifetimes, the rate should closely follow 
the star formation rate (SFR), offering a direct and in- 
dependent way of measuring the SFR and at the same 
time the metal enrichment rate. While this method is 
straightforward in principle, it is hampered particularly 
by the need for corrections due to dust extinction in order 
to derive accurate rates. For Type la SNe, determining 
the rate offers a way to estimate the delay time between 
the formation of the progenitor star and its explosion as 
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an SN, which may shed light on the still mostly unknown 
scenarios leading up to the explosion of the degenerate 
white dwarfs that are assumed to be the progenitors of 
Type la SNe (e.g., Madau et al. 1998; Dahlen & Fransson 
1999; Gal- Yam & Maoz 2004; Maoz & Gal- Yam 2004). 

Due to their potential as standard candles used to mea- 
sure the expansion rate of the universe (Riess et al. 1998; 
Pcrlmutter et al. 1999; Astier et al. 2006), a large num- 
ber of surveys aiming at detecting Type la SNe have been 
conducted during the last decade, resulting in a multi- 
tude of rate estimates. For a recent compilation of rates 
from the literature, see Graur et al. (2011). 

While these results show a fair consensus that the Type 
la SNR increased by at least a factor of ~3 between red- 
shift z = and z ~ 1, there is still an uncertainty in the 
rates at z £ 1, reflecting the uncertainty in such searches, 
the limited statistics and relatively high systematic errors 
concerning searches at z > 1. Moreover, optical searches 
are also hampered by the redshifting of the SN spectral 
energy distribution (SED), making them drop out of op- 
tical filters at z ^ 1. Deep surveys in the near- infrared 
(IR), such as the Hubble Space Telescope (HST) based 
CANDELS (Grogin et al. 2011; Koekemoer et al. 2011) 
and CLASH (Postman et al. 2012) programs, should 
help alleviate this scarcity of high-redshift Type la SNe. 
Rodney et al. (2012) demonstrate the possibility of de- 
tecting and follow-up of Type la SNe at z > 1.5 using IR 
searches. 

While Type la SNe have been the focus of a large num- 
ber of surveys during the last decade, there has been a 
lack of dedicated surveys aiming at finding CC SNe at 
cosmological distances. Most of the existing surveys suf- 
fer from severe selection biases since the foremost aim 
has been to select clean samples of Type la SNe and little 
follow-up has been devoted to non-Type la SN transients, 
including CC SNe. Furthermore, CC SNe are harder to 
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detect and follow on the light curve because they are 
typically ^2 mag fainter than Type la SNe. Therefore, 
there are significantly fewer CC SNRs reported in the 
literature compared to Type la SNRs. 

The local CC SNR (outer distances between 11 and 
^200 Mpc) has been determined by Cappellaro et al. 
(1999), Smartt et al. (2009), Li et al. (2011), Bot- 
ticella et al. (2012), and Mattila et al. (2012), while 
at cosmological distances rates have been determined at 
z = 0.21 by Botticella et al. (2008), z = 0.26 by Cappel- 
laro et al. (2005; these rates are, however, superseded by 
the rates in Botticella et al. 2008), z = 0.3 by Bazin et 
al. (2009), z = 0.66 by Graur et al. (2011), z = 0.39 and 
z = 0.73 by Dahlen et al. (2004), and most recently at 
z = 0.39 and z = 0.73 by Melinder et al. (2012). 

Most rates are consistent with each other and they in- 
crease with look-back time, reflecting the increased SFR 
between z = and 2 = 1. However, there are some 
significant differences in the measured local rates. Fur- 
thermore, there are also claims by Horiuchi et al. (2011) 
that the SNR in general is a factor of ~2 lower than pre- 
dicted from the SFRs, suggesting that there are possible 
unaccounted systematic effects that may affect the de- 
termined rates. Most importantly, necessary corrections 
to account for the dimming of the SNe by dust in the 
host galaxies are highly uncertain. The effect of dust ex- 
tinction is also expected to increase with redshift since 
the rest frame probed in optical searches approaches the 
UV part of the spectrum where the effect of extinction 
is more severe. There should also be a population of CC 
SNe almost completely hidden from optical searches in 
highly extinguished luminous infrared galaxies (LIRGs) 
and ultra-luminous infrared galaxies (ULIRGs; see, e.g., 
Mannucci et al. 2007; Mattila et al. 2004, 2012). Since 
the fraction of the star formation that occurs in U/LIRGs 
increases significantly with look-back time and actually 
dominates the total star formation at z ~ 1 — 2 (e.g., 
Magnelli et al. 2011), we expect the number of CC SNe 
hidden in these environments to increase with redshift. 
The CC SNe hidden in such galaxies must be accounted 
for when the aim is to derive the complete rates of CC 
SNe. 

Also, having a good knowledge of how the SN popula- 
tion is divided into different subtypes and the luminosity 
function of each subtype is essential for deriving accurate 
rates. In particular, knowing the fraction of faint CC SNe 
(M ^ -15) is important when deriving rates since these 
will remain mostly undetected and must be corrected for 
when estimating rates. 

In this paper, we use the SN sample detected in 
the extended GOODS/PANS SN survey conducted with 
HST/ ACS 2002-2005 to derive the rate of CC SNe to 
redshift z ~1.1. In a previous paper (Dahlen et al. 2004), 
we presented CC SNR to z=0.7 based on a subsam- 
ple from the same search (~36% of the search that was 
available at that time). The present paper accompanies 
Dahlen et al. (2008), where we presented the Type la 
SNR to z ~1.6 based on the full sample. This paper is 
organized as follows. In Section 2, we describe the SN 
search and present the sample. In Section 3, we describe 
how SNRs are calculated. The results are presented in 
Section 4 together with a discussion in Section 5. We 
summarize our results in section 6. Throughout this pa- 
per we assume a cosmology with S1m=0.3, 51a=0.7, and 



h=0.7. Magnitudes are given in the Vega system. 

2. DATA 

2.1. Observations 

The data used in this investigation were obtained dur- 
ing HST Cycles 11-13 (2002-2005) and consist of 1 ref- 
erence and 13 search images in the GOODS-North and 

I reference and 9 search images in the GOODS-South 
fields (Giavalisco et al. 2004; Riess et al. 2004, 2007). 
Both fields were observed during two separate sequences 
with a cadence between epochs of *~45 days within each 
sequence. The spacing in time of the search epochs is il- 
lustrated in Figure [TJ Each epoch consists of ~15 point- 
ings with the F850LP filter using the Advanced Camera 
for Surveys (ACS) on board the HST. The field of view 
of ACS is ~11 arcmin 2 , and therefore each search cov- 
ers ^160 arcmin 2 . With a total of 22 epochs, we have 
searched an area equivalent to 1 deg 2 with a time sepa- 
ration of 45 days on average. 

Besides the search filter, F850LP, we also obtained ad- 
ditional photometry in the F775W filter, and in Cycle 

II also in the F606W filter, which allowed us to obtain 
color information of the SN candidates. For a number of 
SN candidates we also obtained spectroscopic follow-up 
using both HST and ground-based facilities. For more 
information on the observations and reduction techniques 
used here, see Giavalisco et al. (2004) and Strolger et al. 
(2004). 

2.2. The Supernova Sample 

In total, 118 SNe were detected during the full search 
campaign. Based on SN spectra, colors, light curves, SNe 
and host galaxy redshifts (spectroscopic when available, 
otherwise photometric redshifts based on optical+NIR 
bands), we classified the candidates into 56 Type la SNe 
and 62 CC SNe (see Strolger et al. 2004 for details re- 
garding the typing of the SNe). For the CC SNe, which 
are the focus of this investigation, we have spectroscopic 
redshifts of the SN or the host galaxy for 53 SNe, while 
we use photometric redshifts for 6 host galaxies. The re- 
maining three SNe are apparently hostless and for these 
we estimate the redshift based on the SN magnitudes, 
colors, and light curve. The number of CC SNe used to 
calculate the rates is 45, excluding from the full sample 
the SNe that are outside the redshift range investigated, 
0.1 < z < 1.3, are fainter than our adopted cut-off mag- 
nitude m z =25.8 (see Section 3.2), or were detected as 
decliners present already in the reference epoch. 

3. CALCULATING RATES 

The SNR is calculated using Monte Carlo simulations 
as described in detail in Dahlen et al. (2004). In sum- 
mary, the simulations assume a given input rate within 
a particular redshift bin and calculates the number of 
detectable SNe for a search with an observational setup 
identical to the HST search, i.e., field of view, temporal 
spacing between epochs, number of epochs, and detec- 
tion limits. A number of properties characterizing the 
SN population go into the simulations, such as the dis- 
tribution of SN peak brightness, light curve shapes, spec- 
tral evolution, the relative fraction of different subtypes, 
and extinction. Rates are thereafter calculated by ad- 
justing the input rate so that the predicted number of 
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detectable SNe matches the number found in the actual 
search within each redshift bin. Below we give a brief 
summary of the assumptions used in this investigation 
when calculating the rates. 

3.1. CC SN Subtypes 

CC SNe are traditionally divided into different sub- 
types based on spectral features or light curve shapes. 
Since each type has different peak luminosity distribu- 
tions and temporal evolution, it is important to correctly 
characterize each subtype as well as determine the rela- 
tive frequency of each type. 

3.1.1. Peak Magnitudes 

We divide the CC SNe into four main types; Type II 
Plateau (HP), Type II Linear (IIL), Type II with narrow 
emission lines (Iln), and Stripped Envelope SNe (SE), 
where the latter includes types lb, Ic, and lib. Further- 
more, Richardson et al. (2002, 2006) note that the IIL 
and SE types have indications of a bi-modality in the lu- 
minosity function and present peak magnitudes for both 
a normal and a bright subpopulation. Here we adopt 
this subdivision into two different subtypes for IIL and 
SE types. Peak magnitudes and dispersion for the six 
types used here are taken from Richardson et al. (2002, 
2006) and are given in Table [TJ 

3.1.2. Subtype fractions 

To estimate the relative intrinsic fractions of the dif- 
ferent subtypes we use the results from Smartt et al. 
(2009) and Li et al. (2011). Smartt et al. present the 
relative fractions in a volume-limited sample (<28 Mpc), 
in total consisting of 92 CC SNe, while Li et al. present 
a volume-limited sample (<60 Mpc) consisting of 133 
CC SNe. Here we combine these results to improve the 
statistics when deriving the fractions. Note that there 
is an overlap of 24 SNe in the two samples, and we take 
this into account when combining the two. The results 
are presented in Table [TJ The IIP type is most com- 
mon (54.8%) followed by the SE SNe (34.0%). Only mi- 
nor fractions belong to the IIL type (6.1%) and Iln type 
(5.1%). However, since that latter two types are more 
luminous on average than the former two, particularly 
Iln, the observed fraction of SNe in a magnitude-limited 
search will likely have higher fractions of these subtypes. 

To estimate the subfraction within the IIL and SE 
types that belong to the bright population, we use the 
number of bright and normal SNe in each of the samples 
used in Richardson et al. (2002, 2006). However, since 
these are basically magnitude-limited samples, we can- 
not use the observed fraction of bright SNe since they 
will be over-represented due to their brightness, which 
could bias the calculation of relative fractions. For the 
SE type, the sample in Richardson et al. (2006) consists 
of 4 bright and 23 normal SNe. Inspecting the distance 
moduli for these SNe reveals that three of the four bright 
SE SNe are at larger distances compared to the normal 
type, while one of them is at a distance where normal 
SNe are also detected. This suggests that 1/24, corre- 
sponding to ~4% of the SE SNe, belong to the bright 
group, although with high uncertainty. For the IIL sub- 
type, 4 objects out of a total of 16 belong to the bright 
subsample in Richardson et al. (2002). Inspecting the 



distance moduli for these objects does not show any sig- 
nificant difference between the bright and normal sub- 
samples. We therefore assume that 4/16, or 25% of the 
IIL SNe, belong to the brighter population. In the fifth 
column of Table [TJ we also include this subdivision of the 
SE and IIL types. We again note that these fractions are 
uncertain. Furthermore, while these are calculated using 
a local sample, it is possible that the relative fraction of 
different types evolves with look-back time due to, e.g., a 
change in metallicity with time. This adds to the uncer- 
tainty. In Section 4 where we discuss systematic errors, 
we investigate how the uncertainty in the fractions may 
affect the derived rates. 

3.1.3. Fraction of faint M > -15 SNe 

Faint CC SNe with magnitudes M > — 15 may be too 
faint to be detected even in local surveys such as the Lick 
Observatory Supernova Search (LOSS; Li et al. 2011), 
which goes out to ~ 200 Mpc. However, if the intrinsic 
fraction of faint SNe could be determined, then it would 
be possible to account for this population when deriv- 
ing the rates. Richardson et al. (2002) argue that the 
fraction of faint Mb > — 15 CC SNe should be at least 
20%, perhaps much higher. This is consistent with the 
results of Cappellaro et al. (1997), who estimate that 
faint 1987A-like SNe should be 10%-30% of the Type 
II SNe. As a conservative limit on the fraction of faint 
SNe, we adopt a value of 20%, consistent with the esti- 
mate of Richardson et al. (2002). The values listed in 
Table [TJ (Columns 2-4) already result in -4% SNe with 
Mb > — 15 when taking into account the dispersions in 
peak magnitudes. To account for a larger fraction of 
faint CC SNe, we add a population of SNe with peak 
magnitude M B = -14.4, similar to the faint SN 1987A. 
In order to make the total fraction of faint SNe 20%, we 
assume that 19% of the SN population is of the faint SN 
1987A type. These fractions are shown in Column 6 of 
Tabic [TJ 

3.1.4. Light Curve and SED 

Each of the four main CC SN types is represented by 
a light curve and a set of SEDs that changes with time 
on the light curve. We base our light curves and spectral 
library on the templates publicly available at the Web 
site of P. Nugent 

3.2. Detection efficiency 

The detection efficiency of the H ST SN search is de- 
scribed in Strolger et al. (2004) and Dahlen et al. (2008). 
In short, two methods were used to derive the detection 
efficiency First, fake SNe of different magnitudes were 
added to the actual images during the real-time man- 
ual searches. The number of added fake SNe was pur- 
posely limited to minimize fatigue by the people doing 
the searches. Second, we used a Monte Carlo method 
to plant thousands of false point-spread functions in the 
data and thereafter tried to recover them using an au- 
tomated detection algorithm. The Monte Carlo test can 
probe a significantly larger parameter space (combina- 
tions of magnitude and host environments) and is there- 
fore complimentary to the real-time test. The Monte 

6 http:/ /supernova. lbl.gov/~nugent/nugent_templates. html 
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Carlo method makes assumptions on the limitations of 
the manual searches, i.e., setting a threshold for detec- 
tion that matches the human searches. It is therefore 
not surprising that the two methods produce similar re- 
sults in detection efficiency. The results from these two 
methods were averaged to give a final detection efficiency. 
The solid line in Figure [5] shows the derived detection ef- 
ficiency. The vertical dashed line shows the cutoff mag- 
nitude used when deriving rates. By applying a cutoff we 
decrease the need for corrections due to incompleteness 
at the faintest magnitudes. Even though this cutoff ex- 
cludes seven SNe, the increase in error due to the lower 
statistics is less than the uncertainties introduced by the 
incompleteness corrections. Also, these are the faintest 
SNe in our sample and may therefore have larger uncer- 
tainties in their type determination. 

3.3. Dust extinction 

Extinction by dust in their host galaxies makes SNe 
fainter and a number of SNe which without extinction 
would be readily detectable may fall below the detection 
limit of the survey. Even if the peak magnitude is above 
the detection limit, a fainter SN spends a shorter time 
on the light curve above the limit and is therefore easier 
to miss during a search. Correcting for the SNe missed 
due to these extinction effects is essential for deriving 
accurate rates. The model used in this investigation to 
calculate corrections for dust extinction in normal spiral 
galaxies is based on the models in Riello & Patat (2005). 
The model assumes different exponential dust and SN 
distributions in galaxies. The effects of the spiral arm 
structure are considered to be negligible and the SNe 
are assumed to be found in host galaxies with random 
inclinations. Furthermore, the model requires the speci- 
fication of a bulge-to-total (B/T) ratio. Here we use B/T 
= 0. The resulting distribution of extinction values can 
be normalized through the face-on central optical depth 
through the model disk galaxy, Ty(0). Mattila et al. 
(2012) show that ry(0)=2.5 gives a good fit to observed 
values of host galaxy extinctions within 12 Mpc. To fur- 
ther check this, we add the sample in Richardson et al. 
(2006) to the observed sample of host galaxy extinctions 
in Mattila et al. (2012). We find that when using the 
same criteria as Mattila et al. (including the observed 
extinctions in galaxies with an inclination less than 60°, 
and excluding two observed SNe with high extinction, 
A v > 3) the mean extinction is (A v ) = 0.42±0.09 for 
our sample of 34 SNe. This is consistent with the pre- 
dicted value (Ay) — 0.44 from the simulations, after us- 
ing the same inclination restriction and excluding a few 
objects (~1%) with Ay > 3 in the high extinction tail of 
the simulated distribution. 

At low redshift, most star formation, and therefore 
also CC SN production, occurs in galaxies with mod- 
erate amounts of dust extinction which we assume can 
be mostly corrected for by the models described above. 
We note, however, that while these models do predict a 
fairly significant tail of high extinction values in inclined 
host galaxies, they only result in about ~ 1% with high 
extinction values Ay > 3 in galaxies with an inclination 
less than 60°. Based on observed data, Mattila et al. 
(2012) estimate that a fraction of about ~15% may have 
such high extinction even in these low to intermediate 
inclination galaxies and are therefore more likely to be 



missed in searches. 

Furthermore, the fraction of star formation taking 
place in highly dust-enshrouded environments, such as 
LIRGs and ULIRGs, increases rapidly with redshift (e.g., 
Perez-Gonzalez et al. 2005; Le Floc'h et al. 2005; Caputi 
et al. 2007; Magnelli et al. 2009, 2011). For example, 
the results of Magnelli et al. (2011) show that at low red- 
shift ;$ 10% of the star formation occurs in U/LIRGs. 
This increases to ~50% at z = 1. A majority of the 
SNe in these environments will suffer from high extinc- 
tion and should be mostly invisible in optical passbands 
(Mannucci et al. 2007; Mattila et al. 2012). 

To account for the SNe that are lost by optical searches 
due to heavy extinction in dusty starbursts, we use the 
results from Mattila et al. (2012), who calculate the 
missing fraction using an approach somewhat similar 
to Mannucci et al. (2007). First, the fraction of CC 
SNe that is hidden in highly dust-enshrouded environ- 
ments within normal galaxies, LIRGs, and ULIRGs is 
estimated. Thereafter the fraction of the total star for- 
mation occurring in these different galaxy types and their 
relative contributions as a function of redshift are taken 
from Magnelli et al. (2011) to estimate the missing frac- 
tion of CC SNe as a function of look-back time. Using 
the results in Mattila et al. (2012), we find that the 
missing fraction increases from ~19% at low redshift to 
~38% at z ^ 1. To take into account the redshift depen- 
dence of the missing fraction, in each of our redshift bins 
we calculate the volume-averaged missing fraction using 
the results in Mattila et al. (2012). In order to correct 
the rates for the missing fraction, results should be mul- 
tiplied by a de-bias factor that can be derived from the 
missing fraction according to the de-bias factor=l/(l- 
missing fraction). 

4. RESULTS 

The CC SNRs are calculated in three redshift bins 
0.1 < z < 0.5, 0.5 < z < 0.9, and 0.9 < z < 1.3, 
centered at z ~0.39, z ~0.73, and z ~1.11, respectively. 
For these bins, we find rates of 3.00tJ;^, 7.39±^;||, and 

9 - 57 -2.80 in units yr^Mpc^lO- 4 ^. The errors are 
statistical only, and systematic errors are discussed in 
the next subsection. The results are also presented in 
Table [H and are plotted in Figure [3] To emphasize how 
the dust corrections affect the results, we also give the 
raw rates without any corrections, the rates after apply- 
ing extinction correction in normal galaxies, and finally 
also the rate after de-biasing for the missing fraction in 
high-dust environments. It is clear that the importance 
of the corrections due to extinction effects increases with 
redshift. In the lowest redshift bin, the final rates are a 
factor of 1.7 higher than the raw rates. This increases to 
factors 2.4 and 2.7 in the intermediate- and high-redshift 
bins, respectively. 

The rates calculated for 8 of the 22 search epochs 
used here were reported in Dahlen et al. (2004). The 
main difference between the investigations is the bet- 
ter statistics in the current results and the implemen- 
tation of the de-bias factor for missing SNe. For the 
0.1 < z < 0.5 bin we found in the earlier eight epoch 
survey a rate 2.511q75, which is fully consistent with the 
new results within statistical errors, both with and with- 
out de-biasing. For the 0.5 < z < 0.9 bin, we found a 
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rate of 3.96±^;g| yr- 1 Mpc _3 10~ 4 ^ in the earlier search. 
These measurements are consistent within la statistical 
errors with the new measurements before applying the 
de-bias factor. After including the de-bias correction, the 
measurements are just outside the la statistical errors, 
leaving the rates in the 2004 paper which were not cor- 
rected statistically lower than the new rates. However, 
adding systematic errors in quadrature to the statistical 
errors makes the two measurement in this bin consis- 
tent within the la errors. For the highest redshift bin, 
0.9 < z < 1.3, there were not sufficient data on the 2004 
sample to derive a rate. This comparison highlights the 
importance of systematic errors when deriving rates since 
these can give a significant contribution to the total er- 
ror budget of the derived rates. In the following, we look 
more closely at the various systematic errors affecting the 
derived SNRs. 

4.1. Systematic Errors 

There are a number of possible sources that may in- 
troduce systematic errors in SNR estimates. Here we 
estimate how different effects may affect the derived val- 
ues. The results are summarized in Table [3l 

Subtype fractions — Since the SN subtypes have different 
peak magnitudes, dispersion, light curves, and SEDs, the 
division into subtypes will affect the "control time" and 
therefore the derived rates. In this paper, we use a divi- 
sion into subtypes based on a combination of the results 
from Smartt et al. (2009) and Li et al. (2011). To 
test the sensitivity to the subtype division, we recalcu- 
late the rates using the Smartt et al. and Li et al. results 
separately. We find that this changes rates by approxi- 
mately ~1%, ~3%, and ~13% in the three redshift bins, 
respectively. In addition, to account for the uncertainty 
in the fraction of bright IIL and SE SNe, we also calcu- 
late rates after assuming that these fractions vary with 
±50%. This introduces a change in the rates by ~1%, 
~2%, and ~7%, respectively. In Table El we have added 
these two sources of uncertainty. 

Faint M > —15 fraction — Based on the discussion in 
Section 3.1.3, we assume that 20% of the CC SNe are 
faint with M > —15. To investigate the dependence on 
the faint fraction, we also derive rates after assuming that 
10% and 30% of the SNe belong to the faint population, 
respectively. This changes the rates by 10%-19%. 

Peak magnitudes — We take into account the dispersion 
on peak magnitudes as given in Table [T] when deriv- 
ing rates. Nevertheless, there are also uncertainties in 
the peak magnitudes which may affect the derived rates. 
Here we use the estimated peak magnitude uncertainties 
in Richardson et al. (2002, 2006) in MC simulations to 
quantify the systematic errors. We find that rates are 
affected by 4%-10% due to the uncertainty in peak mag- 
nitude values. 

Redshift uncertainty — Wrongly determined redshifts for 
the SN may shift objects in or out of redshift bins and 
therefore cause systematic errors in the derived rates. 
Fortunately, we have spectroscopically determined red- 
shifts of the SN or the host for 39 of the total 45 SNe 
in our sample used to derive rates. We therefore do not 
expect that this effect should be large. Nevertheless, we 



use MC simulations to redistribute the redshifts for the 
objects without spectroscopic redshifts by using an error 
distribution based on the photometric redshift fitting. 
We find that the relative errors due to the redshift un- 
certainty is only a few percent, significantly less than the 
statistical errors. Details are given in Table |3] 

Type determination — For our full sample of 118 SNe, 
we have spectroscopic confirmation of the type for 2 
CC SNe and 30 Type la SNe. For the remaining, their 
available photometry in the z, i, and V filters (includ- 
ing non-detections) was fit to model panchromatic light 
curves (made from low redshift analogs), redshifted and 
/c-corrected to the observed passbands. The peak mag- 
nitudes and host extinction were left as free parame- 
ters in the fit, although constrained by the adopted sub- 
type fractions (Section 3.1.2) and extinction distributions 
(Section 3.3). The quality of the light curve fit dictates 
the certainty in the SN type, which we broadly quantify 
as bronze, silver, or gold (in increasing certainty), follow- 
ing a prescription described in Strolgcr et al. (2004). 

Melinder et al. (2011) use MC simulations to esti- 
mate the fraction of misclassified SNe for a sample with 
the requirement that the objects are detected in at least 
two epochs and that color information exists for both of 
these epochs. They find that 5%-10% of the SNe may be 
misclassified. In this investigation, we similarly require 
that the objects are detected during at least two epochs 
and that color information is available for both epochs. 
While the Melinder et al. (2012) investigation uses R and 
I as detection bands and we use % and z bands, we note 
that both investigations probe similar rest-frame bands 
(B and V) at the mean redshifts for the two surveys, 
respectively. Based on this, it is reasonable that the mis- 
classification rate should be similar for the two surveys. 
However, while the simulations assume that only pho- 
tometric redshift information is available, we have spec- 
troscopic redshifts for a majority of the objects or hosts, 
which should improve type determination. Furthermore, 
for a fraction of our sample, we have an additional color 
available, which also should improve typing. Therefore, 
our misclassification rate should be on the lower side of 
the 5%-10% range. To be conservative, we adopt a mis- 
classification fraction that is in the middle of the esti- 
mated range, i.e., 7.5%. This number agrees with the 
fraction of CC SNe misclassified as Type la SNe in the 
full redshift z < 1.5 sample in Melinder et al. Assuming 
that 7.5% could have a wrong type determination, we 
run MC simulations to investigate how this affects the 
derived rates. We find that rates are affected by 6%- 
13% depending on the redshift bin (see Table [3]). This 
should not dominate over statistical errors. 

Extinction corrections — We use results based on the 
models of Riello & Patat (2005) to derive the host galaxy 
extinction corrections. Based on results in Mattila et 
al. (2012), we normalize the model galaxy extinction to 
Ty(0) = 2.5. In Section 3.3 we showed that this results 
in a mean absorption Ay = 0.44 that is consistent with 
the mean extinction in a sample of 34 observed CC SNe, 
(Av)=0.42±0.09. To account for the uncertainty of the 
observed mean, we renormalize the extinction distribu- 
tions so that the resulting extinction distribution fits the 
mean ±1<7/4 V \ and thereafter rederive the SNRs after ap- 
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plying the new extinctions. We find that these changes in 
the host galaxy extinction distributions modify the rates 
by 5%-15%, depending on the redshift bin. Details are 
listed in Table 02 

Extinction laws — To calculate the wavelength depen- 
dence of the extinction, we use the prescription from 
Calzetti et al. (2000) with i?y=4.05. To investigate how 
sensitive results are to the specifics of the chosen extinc- 
tion law, we also derive results after changing the ratio 
of total to selective extinction, Ry by ARy = ±2. We 
have also tried a Cardelli et al. (1989) law with i?y=3.1. 
Using these different scenarios, we find relative changes 
less than 7%, 9%, and 12% in the three redshift bins, 
respectively. 

De-bias correction — In Section 3.3 we discuss the frac- 
tion of CC SNe missed in highly obscured environments, 
particularly in LIRGs and ULIGSs, and the de-bias 
factor that has to be applied to correct for the missing 
fraction of SNe. Using the uncertainty in the missing 
fraction given in Mattila et al. (2012), we find that 
the resulting rates change by 10%-42%, depending on 
redshift. 

After summing the different contributions to the sys- 
tematic errors we find that statistical errors dominate 
in the lowest redshift bin, while the systematic errors 
are similar to or larger than the statistical errors in the 
two higher redshift bins (see Table 02)- This emphasizes 
the importance of understanding the sources of possible 
systematic errors. In particular, corrections due to dust 
extinction and the fraction of faint M > —15 SNe need 
to be better understood in order to better control sys- 
tematic errors. 

5. DISCUSSION 

5.1. Comparison to rates from the literature 

A compilation of rates from the literature is plotted in 
Figure BJ where we also show the SFR estimates from 
Horiuchi et al. (2011) and Magnelli et al. (2009). The 
latter has been put on the same scale by assuming a con- 
version factor between star formation and CC SNR that 
we calculate in Section 5.2. The figure clearly shows an 
increase in the CC SNR with look-back time, consistent 
with expectations from the evolution of the cosmic SFR. 
However, there is a fairly large discrepancy among the 
five presented local or low redshift rates. The highest lo- 
cal rate is presented in Botticella et al. (2012), followed 
by Mattila et al. (2012), Smartt et al. (2009), Li et al. 
(2011), and Cappellaro et al. (1999). We first note that 
the three first rates are based on the number of detected 
SNe within a specified radius during a temporal base 
line. These rates do not apply any corrections for SNe 
that could be missed due to high dust extinction, low 
intrinsic absolute magnitude of the SNe, or SNe missed 
due to insufficient coverage of the sky (e.g., SN exploding 
in areas that could not be observed due to solar avoid- 
ance or because SNe could also remain undetected due to 
the searches not covering their host galaxies such as low- 
luminosity dwarf galaxies). The quoted rates in these 
surveys are therefore lower limits. The rates of Mattila 
et al. (2012) and Botticella et al. (2012) are consistent 
within statistical error bars. This is not unexpected since 



the two investigations use similar selections: Botticella et 
al. use 14 SNe detected within 11 Mpc during 13 years 
(1998-2010), while Mattila et al. use 24 SNe detected in 
the range 6-15 Mpc during 12 years (2000-2011). The 
latter investigation excludes SNe within 6 Mpc so that it 
would not be biased by the local overdensity in star for- 
mation, and therefore also in SNR. The rate in Smartt 
et al. (2009) is based on 92 SNe within 28 Mpc detected 
during 10.5 years (1998-2008), i.e., it covers a volume 
of ^16 and ^7 times larger compared to the searches in 
Botticella et al. (2012) and Mattila et al. (2012). Simi- 
lar to those investigations, the rate presented by Smartt 
et al. (2009) is not corrected for any detection biases 
and should therefore also represent a lower limit. When 
comparing to the more nearby rates, the Smartt et al. 
(2009) rate is about a factor of ~1.7 lower. This could 
be partly due to the local over density affecting the most 
nearby rates. Also, the larger radius in the Smartt et al. 
survey can contribute to this decrease due to increased 
distance modulus, e.g., an SN at r = 28 Mpc is ~1.4 
mag fainter than at r = 15 Mpc and ~2.0 mag fainter 
than at r — 12 Mpc, this could make intrinsically faint 
or dust extinguished SN drop below detection limits at 
larger radii. The sample of Mattila et al. (2012) be- 
tween 6 and 15 Mpc includes at least two SNe with high 
enough extinction that they would likely have remained 
undetected at a larger distance. Also, even though the 
local universe has been well monitored during the last 
decade, possible systematic effects due to non-coverage 
of different regions of the local volume should increase 
with radius. 

At even larger distances, Li et al. (2011) use a sample 
of 440 CC SNe detected in the LOSS. The SNe in this 
sample, which is the largest of any samples used to de- 
rive rates, are detected to distances of ^200 Mpc. While 
the previously described rates consist of compilations of 
SNe detected over the full sky from various sources, the 
LOSS is a dedicated survey where rates are derived from 
imaging a sample of ^10,000 galaxies during a period of 
slightly more than 10 years. The rates are first calculated 
normalized to the surveyed galaxies if -band luminosity 
in units SNuK (where 1 SNuK corresponds to one SN per 
100 yr per 10 10 i Q (K)). To decrease the selection effects 
against SNe in highly inclined hosts, galaxies with an 
inclination >75° are not included when calculating the 
normalized rates. After deriving the normalized rate, the 
volumetric rate is derived by multiplying by the if -band 
luminosity density derived independently by Kochanek 
et al. (2001). When doing this, Li et al. take into ac- 
count that both the normalized rates and the luminosity 
function differ between early-type and late-type galax- 
ies. Furthermore, they also take into account the rate- 
size relation, which suggests that the normalized rates 
are related to the luminosity or the mass of the host 
galaxy. While the rate derived by Li et al. (2011) is out- 
side the statistical errors of the rates derived at smaller 
distances, we do note that the risk of missing SNe due 
to intrinsic faintness or dust extinction, should increase 
with distance, possibly biasing rates to lower numbers. 
Also, Li et al. point out that the lack of corrections due 
to inclination effects (except excluding galaxies with in- 
clination >75°) leads to the largest uncertainty in their 
rates, possibly underestimating the results. There is also 
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an additional possible systematic effect in that the rate 
calculation requires the use of the local -fT-band luminos- 
ity density. 

The last local rate is from Cappellaro et al. (1999) 
and is built on a compilation of different surveys, many 
of which are based on photographic plates, but also the 
visual search by Robert Evans (Evans 1997). The rates 
are calculated normalized to the B-band luminosity and 
volumetric rates are derived by multiplying by the lo- 
cal -B-band luminosity density (Cappellaro et al. 2005), 
different from Li et al. (2011), who use the if-band. 
Considering the different search techniques, it is almost 
surprising that the derived rate is so close to the rates 
derived from CCDs. However, it is hard to evaluate dif- 
ferent possible systematic effects which may affect the 
derived rates. This is outside the scope of this paper, 
but we do leave the rate in this paper because it is one of 
the first based on a compilation consisting of a significant 
sample of SNe. 

When looking at the local rates, it is obvious that there 
is a trend of decreasing rates as the samples are selected 
at increasing radii, with deviations between some of the 
rates being larger than statistical errors. The trend sug- 
gests that more SNe may be missed at larger radii, pos- 
sibly due to extinction in the host galaxies, which, if un- 
corrected for, may underestimate the rate. There could 
also be a relatively high fraction of intrinsically faint SNe 
that would be more easily missed at large radii, which, 
if not accounted for, would underestimate the derived 
rates. 

One should also note that there is a local overden- 
sity in star formation in the very nearby universe, which 
could drive the most local rates to high values. Assuming 
that the Botticella et al. (2012) rate, measured within 
11 Mpc, is affected by the local over density with ap- 
proximately a factor of 1.5, which is the ratio of the 
nearby SFR within 11 Mpc from 11HUGS (Kennicutt 
et al. 2008; Bothwell et al. 2011) and a mean of a sam- 
ple local measurement at larger radii going out to ^ 100 
Mpc reported in Horiuchi et al. (2011), we find that after 
correction the Botticella et al. rate is still larger than the 
rate from Smartt et al. (2009) et al. and Li et al. (2011), 
even though it still represents a lower limit. On the other 
hand, Botticella et al. (2012) note that when normaliz- 
ing their rate to the local 5-band luminosity density, 
they find a value consistent with the normalized rate at 
z ^0.01 in Cappellaro et al. (1999), suggesting that lo- 
cal overdensity can fully explain their elevated local rate. 
Moreover, Mattila et al. (2012) noted that at very small 
radii, r ^ 6 Mpc, there is an elevated SNR, reflecting the 
over abundance in the SFR in the nearby universe. They 
show that the rate calculated within r <6 Mpc is about 
three times the rate at r <10 Mpc. To avoid being af- 
fected by the high number of SNe found at small radii, 
the rate in Mattila et al. is based on SNe found in the 
range 6 Mpc< radius <15 Mpc, which should limit the 
influence of a local overdensity on their rates. Still, it 
cannot be ruled out that the rate in Mattila et al. is 
elevated due to the effects of the local overdensity. 

A non-local rate at z ~ 0.21 based on ~45 CC SNe 
is presented in Botticella et al. (2008). The SNe are 
from the Southern inTermediate Redshift ESO Super- 
nova Search (STRESS) and the rates replace the pre- 



liminary results from a subsample of the survey given in 
Cappellaro et al. (2005). While this is consistent with 
the local rates in Botticella et al. (2012) and Mattila et 
al. (2012), there is an indication that the rate is on the 
low side to be consistent with the expected increase in 
the SFR, and therefore also SNR with look-back time. 
We note that Botticella et al. (2008) correct for extinc- 
tion in normal galaxies, but do not make any correction 
for U/LIRG like environments, which could have some 
effect on the rates. There could also be a number of sys- 
tematic effects affecting the rate. This includes the fact 
that only about one-third of the transients in the survey 
has an SN-type determination and that for another third 
there is an ambiguity whether the transient is an SN or an 
active galactic nucleus. Furthermore, this was a pointed 
survey, observing a sample of galaxies, with the aim at 
deriving the B-band luminosity-normalized SNR. Deriv- 
ing the volumetric rate based on this includes the extra 
systematic uncertainty when multiplying the luminosity- 
normalized rate with the global luminosity density at the 
redshift range of interest. 

At somewhat higher redshift z ~0.3, Bazin et al. 
(2009) present the CC SNR, based on a sample of 117 
SNe from the Supernova Legacy Survey, which is the 
largest non-local sample of SNe used to derive rates. Sim- 
ilar to the Botticella et al. (2008) rate, this rate is consis- 
tent with lower redshift rates; however, it does not show 
an increase compared to the most recently published lo- 
cal rates, which would be expected if the rate increases 
with look-back time. Again, we note that while the Bazin 
et al. (2009) rate includes extinction corrections in nor- 
mal host galaxies, it does not correct for highly dust ex- 
tinguished regions, which could somewhat bias the rates 
to lower values. Furthermore, there are systematic effects 
that are difficult to quantify in the derived CC SNR. In 
particular, Bazin et al. (2009) do not directly derive 
the CC SNR; instead they estimate the relative rate be- 
tween CC SNe and Type la SNe at z ^0.3 and then 
use the Type la SNR derived at z ^0.5 by Neill et al. 
(2006), which is extrapolated to z ~0.3, and thereafter 
multiplied by the assumed ratio to derive the CC SNR. 

At higher redshifts, overlapping the range presented in 
this paper, Graur et al. (2011) derive the rate at z ^0.66 
based on three SNe classified as CC SNe, while Melinder 
et al. (2012) present rates at z ~0.39 and z ^0.73 us- 
ing three and five CC SNe in the two bins, respectively. 
Both of these results apply corrections for host galaxy 
dust extinction as well as galaxies missing in U/LIRGs, 
using methods similar to those used here. As can be seen 
in Figure 0] these rates are consistent with the rates de- 
rived here in the overlapping redshift range. At the same 
time, the error bars on the Graur et al. and Melinder et 
al. rates are large, covering a factor 2-10 in the derived 
rates, mainly due to the low SN statistics. Additional 
systematic errors could also be large in those results due 
to, e.g., the lack of spectroscopic confirmation of the red- 
shifts of both the SNe and the host galaxies. With our 
new rates, based on the largest high-redshift statistical 
sample for which we have spectroscopic confirmation of 
the redshift of the SN or host for a majority of the sam- 
ple, we are able to significantly decrease the size of the 
error bars, making it possible to set firmer constraints on 
the evolution of the cosmic SNR. 
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5.2. The star formation rate connection 

One of the most important diagnostics that the CC 
SNe can be used for is the measurement of the cosmic 
SFR. Note that the CC SNR provides a completely inde- 
pendent method for measuring the evolution of the SFR, 
compared to other widely used methods based on, e.g., 
galaxy luminosity densities (UV, IR, Ha, etc.). For a di- 
rect comparison between the local CC SNR and different 
SFR indicators, see Botticella et al. (2012). Previous 
results have shown that the increase in CC SNR with 
look-back time is consistent with an increase in the rate 
as expected from the SFR (Dahlen et el. 2004). How- 
ever, the details in the evolution of the CC SNR and their 
normalization to the SFR still have large uncertainties, 
which has prevented a more detailed comparison. 

In Figure HI we plot the CC SNR with the SFR from 
Magnelli et al. (2009) and Horiuchi et al. (2011). To 
relate the CC SNR and the SFR, we assume a Salpeter 
initial mass function (IMF; Salpeter 1955), tp(M), in the 
range 0.1 < M/Mq < 125 together with the assumption 
that all stars with masses in the range 8 < M/Mq < 
50 explode as CC SNe. We can thereafter relate the 
SFR (SFR in units of Moyr^Mpc' 3 ) and the CC SNR 
(SNR in units of yr" 1 Mpc~ 3 ) by 

SNR(z) = kxh 2 x SFR(z), (1) 

where 

Q(M) rf Af 

ioi 5 M: M ^ dM 

Numerically, this results in fc=0.0070 MX . There are 
two main uncertainties in deriving this value. First, the 
choice of IMF affects the result. However, the depen- 
dency of the IMF in Equation (2) is mostly cancelled out 
since the derived SFR in Equation (1) also depends on 
the IMF. E.g., Melinder et al. (2012) show that even if 
making a significant change in the IMF, such as choos- 
ing a modified Salpeter IMF from Baldry & Glazebrook 
(2003), which will increase k by a factor of ~2, the re- 
sulting change in the relation between SNR and SFR is 
only ~2% due to the change of the SFR normalization. 
The second uncertainty comes from the selected range of 
CC SN progenitor masses. Here we use a lower limit of 
8Mq consistent with the results of Smartt et al. (2009), 
who estimate that the minimum stellar mass for type 
IIP progenitors lies in the range 7 < M/Mq < 9.5. This 
range corresponds to an uncertainty in k by +22%/-23%. 
The upper limit (50M Q ) is taken from Tsujimoto et al. 
(1997), who suggest that progenitors more massive than 
this become black holes directly instead of exploding as 
SNe. Due to the steep slope of the IMF, the actual choice 
of upper limit has limited impact on the derived fc value. 
Here we assume that the upper limit is restricted to the 
range 30 < M/Mq < 100, resulting in an uncertainty 
of +6%/-9%. Taking these ranges in upper and lower 
limits into account, we find that k should lie in the range 
0.0048< k <0.0089, indicating a maximum uncertainly 
in the normalization by less than ±25%. We also note 
that there is an uncertainty in the normalization of the 
SFR by approximately ~20% (e.g., Figure 5 in Horiuchi 
et al. 2011). An uncertainly of this order is also illus- 
trated on Figure [3] as the difference in SFR between the 



estimates in Magnelli et al. (2009) and Horiuchi et al. 
(2011). Therefore, there should still be a total systematic 
uncertainty of ~32% when comparing the SFR derived 
from galaxy luminosities and the SFR from CC SNe. 

Using the relation in Equation (1), we see in Figure |4] 
that our newly derived rates well match the predictions 
from the SFR over the range 0.4 ;$ z ^ 1.1. Indicating 
that the assumptions made when deriving the rates are 
not severely affected by systematic errors, and the as- 
sumed relation between SFR and SNR seems reasonable. 
We note that the local rates by Botticella et al. (2012) 
and Mattila et al. (2012) are also consistent with the 
expectations from the SFR, even though these rates rep- 
resent lower limits. 

Based on the CC SNRs published prior to 2011, Hori- 
uchi et al. (2011) noted that the observed CC SNRs and 
the rates expected from the SFR are systematically off- 
set by a factor of —2, particularly at lower redshifts. The 
main explanations for this discrepancy suggested that 
a high fraction of CC SNe could be missing from the 
searches due to a significant population of intrinsically 
faint SNe or that a large fraction is invisible to searches 
since they are hidden behind significant amounts of dust. 
In this investigation, we take both of these sources into 
account when deriving rates. Since our results agree well 
with what is expected from the SFR, it seems likely that 
the discrepancy in Horiuchi et al. (2011) is indeed due to 
the factors suggested by the authors. For the other high- 
redshift rates, the inclusion of corrections for SNe missed 
in U/LIRGs in the rates of Graur et al. (2011) and Melin- 
der et al. (2012) should contribute to increase these rates 
to make them consistent with the SFR. For local rates, 
Mattila et al. (2012) report a value that is also consistent 
with the SFR. They argue that within 15 Mpc, a large 
fraction of either faint or dust-enshrouded SNe should 
not be missing, which is the reason that their rate is 
consistent with the SFR. Similar arguments should hold 
for the rate in Botticella et al. (2012), even though this 
rate could be elevated more due to the local overdensity 
in star formation. For the remaining rates in Figure |4] 
that fall below the SFR lines, it is likely that these two 
effects both contribute to the rates, in combination with 
unaccounted systematic errors that could affect each dif- 
ferent survey. 

We finally note that the good agreement between the 
SFR and the CC SNR derived here over the whole red- 
shift range to z —1.1 indicates that we are not signif- 
icantly affected by incompleteness at higher redshifts 
where the SNe become fainter. Specifically, this demon- 
strates that we are quite complete to m(f8501p)=25.7, 
typical of a CC SN at z -1.1 (with M B =-17.4). Sim- 
ilarly, for the Type la SNR published in Dahlen et al. 
(2008), based on the same SN search, we should there- 
fore also be complete to m(f8501p)=25.7, corresponding 
to faint SNe la (with M B =-18.3) at redshifts z ~1.6, the 
highest bin where the Type la rate was derived. This sug- 
gests that the drop in the Type la rate at z > 1.4 noted 
in Dahlen et al. (2008) is not due to incompleteness. 

6. CONCLUSIONS AND SUMMARY 

We have used a sample of 45 CC SNe, out of a total 
of 62 detected with ACS during the #ST-extended SN 
search in Cycles 11-13, to derive the CC SNR to z ~1.1. 
Our main conclusions are as follows. 
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• After correcting for host galaxy extinction and 
the missing fraction in highly obscured environ- 
ments, we find rates 3.00±J;|f ±J;§y at (z)=0.39, 



discrepancy between some earlier derived 



0.3 



7 oq + 1-86+3.20 „, 
' - Ja -1. 52-1.60 dl 



(z}=0.73, and 9.57+£ 



+3.76+4.' 



80-2.80 



at 



(z) =1.11. The rates are given in units of 
yr _1 Mpc~ 3 10~ 4 hj Q . The first errors represent 
statistical while the second ones are the estimated 
systematic errors. 

• Our rates at z ^ 0.4 are consistent with those ex- 
pected from the cosmic SFR. 

• Statistical uncertainties dominate our rate esti- 
mates in the low-redshift bin, however; systematic 
errors start to dominate at z > 0.5. 

• While the most recent local CC SNRs are consis- 
tent with those expected from the SFR, there is a 



rates and the SFR. As noted in Horiuchi et al. 
(2011) and Mattila et al. (2012), this is likely 
due to a combination of SN missed due to dust 
extinction and/or faint intrinsic luminosity and un- 
accounted systematic effects. 
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FlG. 1. — Cadence of the 13 search epochs in GOODS-North and 9 search epochs in GOODS-South. 
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Fig. 2. — Distribution of detection magnitudes for the full sample of 62 CC SNe found in the HST search. The solid line shows the 
detection efficiency derived by adding fake SNe to the search images. The vertical dashed line shows the magnitude cutoff applied when 
calculating rates. 
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FlG. 3. — Core collapse SNR. The circles show "raw" rates without any correction for extinction or "missing fraction". The triangles 
show rates after correcting for extinction in normal galaxies. The squares show our final rates after also correcting for the fraction of CC 
SNe expected to be missed in dust-enshrouded environments, in particular in LIRGs and ULfRGs. 
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Fig. 4. — Core collapse SNRs from this work together with rates from the literature. Shown also are the star formation rates from 
Magnelli et al. (2009) and Horiuchi et al. (2011) after assuming a relation between star formation and CC SNR given in Section 5.2. Rates 
from the literature are shown with statistical errors. For the rates in this work, we show statistical errors with black error bars, while 
gray (larger) error bars show the added statistical and systematic errors (in quadrature). Note that the local rates within 11-15 Mpc by 
Botticella et al. (2012) and Mattila et al. (2012) may be elevated due to the local overdensity in star formation. See Section 5.1 for details. 
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TABLE 1 

Subtypes of Core Collapse Supernovae 



Type 


Peak M B 


Dispersion 


Fraction™ 


Fraction' 1 


Fraction 


IIP 


-16.67 


1.12 


54.75% 


54.75% 


44.31% 


IIL normal 


-17.23 


0.38 


6.10% 


4.57% 


3.70% 


IIL bright 


-18.94 


0.51 




1.52% 


1.23% 


Iln 


-18.82 


0.92 


5.12% 


5.12% 


4.14% 


SE normal 


-17.07 


0.49 


34.03% 


32.61% 


26.24% 


SE bright 


-19.38 


0.46 




1.42% 


1.31% 


Faint (M B > -15) 


-14.4 


0.5 






19.06% 



Note. — "Fractions of main types derived from data in Smartt et al. (2009) and Li et al. (2011). ^Fractions after dividing the SE and 
IIL subtypes into a normal and a bright population, based on the data in Richardson et al. (2002; 2006). c Fractions after assuming that 
in total 20% of the SN population is fainter than Mb >-15. 



TABLE 2 

Type Core Collapse Supernova Rates 



Redshift Bin 


Redshift" 


SNR 6 


Error(stat+sys) 


SNR C 


Error (stat+sys) 


SNR d 


Error (stat+sys) 


N 


N dist 


O.K z < 0.5 
0.5 < z < 0.9 
0.9 < z < 1.3 


0.39 
0.73 
1.11 


1.79 
3.14 
3.51 


+U.77 +0.30 

-0.56 -0.26 
+0.79 +0.70 
-0.64 -0.52 
+ 1.38 +0.82 
-1.03 -0.91 


2.24 
4.86 
5.95 


+0.U6 +0.46 

-0.70 -0.37 
+ 1.22 +1.36 
-1.00 -0.92 
+2.34 +1.81 
-1.74 -1.60 


3.00 
7.39 
9.57 


+ 1.28 +I.U4 

-0.94 -0.57 
+ 1.86 +3.20 
-1.52 -1.60 
+ 3.76 +4.96 
-2.80 -2.80 


9 
25 
11 


9.99 
23.56 
11.44 



Note. — "The effective redshift is defined as the redshift that divides the volume in the redshift bin into equal halves. 'Rates without 
any corrections. c Rates corrected for extinction in normal galaxies. d Final rates corrected for SNe hidden in high extinction regions, 
particularly in U/LIRGs. Rates are given in units yr _1 Mpc~^ 10~ 4 h'^ , assuming a cosmology with S7jvf = 0.3 and Qa = 0.7. Statistical 
and systematic errors are shown separately. The last two columns give the number of SNe in each bin, where N is the "raw" counts and 
N dist gives the number after taking into account the redshift probability distribution of the SNe. 



TABLE 3 
Error sources 





Redshift 


Redshift 


Redshift 




O.K z <0.5 


0.5< z <0.9 


0.9< z <1.3 


Subtype fraction 
Faint (M > -15) fraction 
Peak magnitudes 
Redshift uncertainty 
Type determination 
Extinction correction 
Extinction laws 
Missing fraction 


+0.9% 
-1.4% 
+ 13.5% 
-10.0% 
+4.7% 
-4.0% 
+4.1% 
-2.7% 
+7.5% 
-9.2% 
+9.7% 
-4.6% 
+6.5% 
-6.3% 
+28.0% 
-9.6% 


+3.7% 

-4.2% 
+ 18.6% 
-13.2% 
+8.7% 
-7.7% 
+ 1.1% 
-1.8% 
+7.8% 
-6.4% 
+ 14.5% 
-6.9% 
+8.8% 
-4.0% 
+33.0% 
-10.3% 


+ 13.9% 
-15.0% 
+ 14.0% 
-13.5% 
+ 10.5% 
-9.7% 
+2.5% 
-1.7% 
+6.1% 
-13.2% 
+ 15.5% 
-6.8% 
+ 11.8% 
-0.2% 
+42.0% 
-11.6% 


Systematic summed 


+34.6% 
-19.0% 


+43.3%, 
-21.6% 


+ 51.8% 

-29.3% 


Statistical errors 


+42.8% 
-31.2% 


+ 25.2% 
-20.5% 


+39.3% 
-29.2% 



Note. - Different sources contributing to systematic uncertainties. For the summed errors, the difference sources are added in 
quadrature. See the text for details. 



